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SUMMARY 


Methods are described which may be used to transform the equations 
of motion of a vehicle to a standard form, measure the state of the system 
when described in this standard form, and remove first flexure mode effects 
from pitch rate signals. The methods are described by the use of illustrative 
examples. 


INTRODUCTION 

In state vector control of linear systems it is assumed that the differential 
equations of motion for the plant can be written as n linear first-order equations 
where n is the order of the system. If time-optimal or average power optimal 
control is to be used, each of the n state variables must actually be measured. 

When the plant is an aerial or space vehicle, the differential equations of 
motion that serve as the starting point in the control synthesis problem are 
neither first order nor written in terms of state variables that can be measured 
with existing sensors. Because of this, it is necessary to apply certain trans- 
formations to the original equations. A procedure will be demonstrated on the 
pitch equations of motion for a flexible aerial vehicle. It will be assumed that 
gyros are to be used as the primary measuring devices. Accelerometers or 
strain gages could be used as well. 

There are three transformations to be considered: 

• A transformation to put the equations of motion in a standard form 
involving no derivatives higher than first 

• A transformation to permit measurement of the body bending dis- 
placements 

• A transformation to permit the measurement of angle of attack with- 
out an external probe 

Then a scheme is described which may be used to obtain a measurement 
of a vehicle pitch rate that is free from the effects of the first flexure mode. 
This scheme involves blending the signals from two rate gyros located fore 
and aft of the first bending antinode of the vehicle. 


The method for measuring the state vector is from an internal report 
of the Minneapolis-Honeywell Regulator Company (ref. 1). The procedure 
for obtaining a pitch rate signal free from first mode flexure effects is taken 
from ref. 2. 


SYMBOLS 


CG Center of gravity 

M Pitching moment (positive up) 

U Forward speed 

X A force along the x axis (positive forward) 

Z A force along the z axis (positive down) 

g Gravitational constant 

q = 9 Pitch rate 

w Velocity along the z axis (positive down) 
x Distance along the longitudinal axis (positive forward) 

z^ j Deflection of the ( ) body bending mode 

^ By. definition 

a Angle of attack 

6e Elevator deflection 

6T Throttle deflection 

9 Pitch angle (positive nose up) 

Subscripts 

6e, 6T, U, a, a, q, w, z . , z^, z., 9 indicate partial derivative. 
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STANDARD FORM FOR THE EQUATIONS OF MOTION 


As stated before, the differential equations describing a linear plant which 
is being forced by r forcing functions can always be written as an equivalent set 
of n linear first-order differential equations of the form 


dx 


l = : 


dt 


= x, = a„,(t)x 1 + a 10 (t)x 0 + . . . + a, (t) x_ +b 11 (t)u. 


1 111 12 ' ' 2 


In w n “11 w “1 


+ . . . + b (t) u 
lr r 


dx 9 

-at = x 2 = a 21 (t) *1 + a 22 (t) *2 + ' • ' + a 2n (t> x n + b 21 (t) U 1 


+ . . . + b 2r <t) u r 


dx 

-3——- = x = a 1 (t) x 1 + a 9 (t) x 9 + . . 
dt n nl 1 n2 2 


• + M X r> + b „1 M U 1 

nn n n 1 1 


+ . . . + b (t) u 

nr r 


( 1 ) 


The x^'s are referred to as the state variables, and the transformation given 
by equation ( 1 ) is called the point transformation. When the system equations 
are written in this form they will be said to be in standard form. This point 
transformation can be written as a vector differential equation of the form 


dx _ 
dt 


Ax + Bu 


( 2 ) 


To demonstrate how a set of equations is put in the form of equation ( 1 ) two 
simple examples are presented. Besides being useful in themselves, the two 
examples should be sufficient to demonstrate the general procedure and provide 
motivation for treatment of the aeroelastic vehicle which is presented later. 
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A Single-Degree-of-Freedom System 

The differential equation of motion for a fourth-order single-degree-of- 
freedom system can be written as 

*x + a 3 'x + a 2 ’x + ajX + a Q x = bu (3 


A 

X 1 - X 

A • 

x 2 X 1 

A • 

X 3 ‘ x 2 

A • 

x 4 ~ x 3 

equation (3) can be written as a system of first-order equations: 


-a -a, 

o 1 


When written as a vector differential equation, equation (5) becomes 


x = Ax + Bu 


where 


” 3, rj ”&r 


k 



Three-Degree-of -Freedom Equations 
for Longitudinal Motion of a Rigid Aircraft 

More effort is generally required to reduce multi -degree-of -freedom 
systems to the standard form than is required for single-degree-of -freedom 
systems. The pitch equations of motion for a rigid aircraft illustrate the 
difficulty without requiring matrix manipulations. 

The equations of motion are 


U - X V V - X o a + g9 


a 


X 6t 6T = 0 


u. 


Z TT U + a - Z a -9 
U w 


V- Z 6e 5e ' 0 
o 


(7) 


- M„U - M.« 


M a + 9 

a 


M 9 

q 


M 6e 6e ' M 6T 6T = 0 


In the single-degree -of -freedom case, the equation could be reduced to 
standard form by simply defining the derivatives as new variables. Trying 
this again, let 


X 1 ‘ 9] x 2 ~ X 3 " a; X 4 “ U ’ U 1 " 6e ’ U 2 ' 6 T 


Then equations (7) become 


' X u x 4 ‘ X a x 3 + ^1 X 6T U 2 = 0 


( 8 ) 


"U Z U X 4 + x 3 Z w X 3 “ X 2 


U 


Z 6e U 1 = 0 


(9) 


^4 - M a k 3 


M <* x 3 + x 2 " M q X 2 ' M 6e U 1 ' %T u 2 


= 0 


The first and second equations are well set; they can be solved for x 4 and x„ 
in terms of zero-order derivatives. The third has two first-derivative terms; 
therefore neither Xo nor x» can be solved for in terms of zero-order deriva- 
tives. Since there is already an equation for x„, the obvious solution is to 
eliminate x„ from the third equation by use of tne second equation. The third 
equation then becomes 
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- M U ' M a U o x 4 + \' M a " M a Z w x 3 + *2 + " M q ' M a x 2 + 


■ M 6e - u l 

o / 


M 6T u 2 = 0 


Now the equations can be written out by using the first of equations (8) the 
first two equations of (9), and equation (10). 


22 *23 a 24 x 2 


l 32 a 33 a 34 x 3 


b 21 b 22 I U 1 


b 31 0 


a 41 0 


a 43 a 44 x 4 


where 


M. + M 
a q 


M. Z + M 
aw a 


M -f M. 

1VA U a U 


b 91 = M» - + Mr 

21 a U Oe 

o 


a 32 " 1 


a 44 = X U 


It should be clear from this example that had all derivatives appeared in every 
equation, it would have been more difficult to put the equations in standard form. 


6 



Elimination of Derivatives of the Forcing Function 


In some instances the right sides of the plant equations naturally contain 
derivatives of the forcing functions. If a multi-degree-of -freedom system is 
written as a single equation by using Laplace transform techniques deriva- 
tives of the forcing function will appear on the right. The right side of the 
differential equation is equivalent to the numerator in Laplace transform no- 
tation. 


Numerator or right-side plant dynamics in conjunction with relay-type 
lumps in the forcing functions would require any analysis or synthesis technique 
to accommodate derivatives of impulse functions This problem can be avoided 
formally by transferring to a set of coordinates in which only zero-order forcing 

terms appear. 


Equation (13) contains derivatives of the forcing function f(t); that is, the 
differential equation is of the form 


,n 

d x 


,n-l m-2 

n + + a 2« t)4 ^f + ' • +a 


dt 


= b <t) 
o 


dt 


dt 


n-l (t) lf + 


a (t) x 
n 


(13) 


d n f 


dt 


n 


+ b i (t ) 




dt 


n-l 


+ b (t) f (t) 
n 


In the case of constant coefficients this equation corresponds to a transfer 
function of the form 


X(s) _ N(s) 

F(sT D(s) 

where N(s) and D(s) are finite polynomials in the Laplace variable s. 

To write equation (13) in standard form free from derivatives of the 
forcing function f(t), the following technique can be used (ref. 3). Let 

x = Xj + G q (t) f 



x 2 + G l (t) f 


(14) 
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dx 


n- 1 


dt 


x 

n 


+ G - (t) f 
n - 1 


dx 

TT" ‘ a n (t), ‘l - a n-l (t,x 2 ' ' a l (t)3t n + G n (t) f 

Now the G.'s can be determined by successively eliminating x., x„, . . . , > 
from (14) and requiring that the resulting differential equationbe Identical 
with (13). This leads to the following recurrence formula for determining 
the G.'s. 6 


G q <*> 

G. (t) 


b o (t) 


i-1 i-k 

b i (t) -t T. 

k= o s= o 


n + s - i| 
n - i 

I 


a. , (t) 

l-k-s 


d s G, 


dt 


(15) 


wherej f is the (n - i + lj th binomial coefficient of order n + s - i. 
These binomial coefficients are defined by the following formula 


(P| = P.' 

[q| jq? tp-q)'.j 

As an example of the use of the recurrence relation (15), consider the system: 

2 • 

e + 2?we+w e = -au - bu 
where a, b, u, are constants. 

Comparing with equation (13), 


n = 3 


a l = 

2?w 

b = 
o 

0 

a 2 

2 

w 

b l " 

0 

a 3 

0 

b 2 

-b 

f(t) = 

u(t) 

II 

CO 

Xi 

-a 
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Equations (14) become 


e = x, + G u 
1 o 


Xi = x 2 + GjU 


*2 ” X 3 + G 2 U 


Xg = -W Xg - 2 ? U Xg + GgU 


According to (15) 


G = b = 0 

o o 


d s G 


'i*> • b i<‘> T <2 S 8 > a i-s-rr = 0 

at 


s= o 
1 2-k 


<*> * b 2 - I I t ”> 


d S G, 


a 


k=o s=o 


2-k-s ,.s 

dt 


b 2 


I <'t S > 


s J) 
dG 'o 


8= O 


2 - s -v 


y (i + s, . 

s 4 0 1 1 a >- 8 ysF 


- b„ = - b 


2 3-k 

°3 <‘> * b 3 - I I <8> a 3- 
k= o s= o 


d S G, 

k-s ~~r = _a ‘ a i g 2 

dt 


= -a + 2 5 u b 
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APPLICATION TO THE PITCH AXIS OF AN AEROELASTIC VEHICLE 


If the pitch axis of an aeroelastic vehicle is to be time-optimally con- 
trolled, it is necessary to be able to measure with standard instrumentation 
the body-bending displacements and the angle of attack. Techniques for doing 
this are presented in the following subsections. 


Measurement of Body-Bending Displacements 

To permit the measurement of the body-bending displacements, with 
standard instrumentation, it is necessary to make three assumptions: 

1. The location of the center of gravity is known. 

2. The mode shapes are known. 

3. Only a finite number, N, of modes exists. 

The first two of these are naturally compatible with a basic assumption of the 
report, i. e. , the characteristics of the vehicle are known. The third assump- 
tion must be made for any analysis or synthesis technique that uses the eigen- 
functions of the partial differential equations that describe the flexure of the 
free-free body under accelerated or unaccelerated motion. Since the eigen- 
functions present the best approximation to the characteristics of a large com- 
plex structure that is likely to be found in the near future, it appears that 
the third assumption is not restrictive. 


Consider figure 1 where 9 G is the pitch attitude sensed at the position x^. 
Then 


0(t) - SqU) + tan * (x^., t) 


(16) 


For 9z/9x small, tan * 9z/9x ~ 9z/9x. Assuming this approximation is 
satisfactory, equation (16) can be rewritten as 


e G (t) = e m - -gj- (x G , t) 


Assumption 3 


is mathematically z(x,t) 



(17) 


Th 



x 


is the normalized z coordinate of the ith flexure mode at station x. 
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Then 9z/3x 


9 g ( ‘> - 


N 



N 

<Mt) - X 

i = 1 


z^(t) so that equation (17) is by assumption 3 



( 18 ) 


Assumption 2 states that(g) x is known and hence also @ x = d/dx @ x 
• _ , 2 N. If (N + 1) altitude sensors are placed at (N + 1) different 

positions aiong the fuselage, an algebraic system °^ ( ^ + 1 h near equaUons 
for the (N + 1) unknowns, 6 , z\, Zo» • • • » Z N» 18 °kt ainec * '*• ®* * 
sensor is placed at x Gk and its output denoted by 0 Gk , the following system 

results): 


e r (t) = e<t) 
g K 



K = 


1 . 2 , 


N + 1 


(19) 


This system may be written in vector form as 

9 q - Az 


( 20 ) 


where 


9 


G 


z 


is a column vector with Kth 
is a column vector with Kth 


component 

component 


GK n 

' 6 for K = 1 ' 

z K - i for K> lj 


and 


A = (a-Jwitha.j 


1 and a i; j + i 



for j >-l 


Note that A is a function of the mode shapes only, 
the sensors can be chosen such that A is nonsingular. 


Thus, locations for 
Then 


z = 



( 21 ) 
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and differentiating with respect to time. 


z 


a-‘* g 


( 22 ) 


Thus if (N + 1) rate sensors are used, d G is obtained, and since A is known 
z may be computed. 


The entire procedure has been developed on the knowledge that each 
sensor in essence measures one unknown and that if sufficient sensors are 
used, sufficient equations can be written to determine the variables of the 
equations of motion. The only problem is to place the sensors so the A 
matrix is nonsingular and well set. This need not be a matter of great con 
cern in that the sensors can be first arbitrarily placed for convenience and 
the resulting A matrix checked for singularity and to see that it is not ill- 
conditioned. 


Measurement of Angle of Attack 

Since the typical angle-of-attack sensor is unacceptable for flight control 
purposes and is likely to remain so for many important applications, it is 
desirable to be able to derive a high«-quality substitute. The angle of attack 
that can be measured without an external probe is that between the relative 
wind and the longitudinal reference line. This proves adequate for flight 
control purposes. 6 


Now consider figure 2 where the x and z axes are shown at time t. 
time vt + A t) they are shown as x' and z 1 axes. 

Let 


At 


P z ~ the z coordinate of P and P 1 = z coordinate of P’ 

z 

Then 

P Z = z(x A' t) 

and 

P z' = ‘Ursine* - <x A - x CG ) sin [ © (t + At ) - 0 (t )j 
^ £z(x a , t t ^^t) cos [ 0<t + At) - 6 <t)]J 

Since 

0<t + At)-e<t) = e<t)At + o[<At) 2 ], 
cos [$(t + At) - e(t)] = 1+0 [(At) 2 ] 
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and 

sin [e(t + At) - 0<t)] = $(t)At + 0 [(At) 2 ], 

then 

— Tt — ~ = ‘ U o sina ^ ’ ^ X A ' X CG* + 

z(x A , t + At) - z(x . , t) 

A — ^ + O(At) 

At 


Therefore, 


w(x A , t) 


- U q sin a (t) - (x A - x qq) 0(t) + 


9z(x a , t) 


at 


( 23 ) 


where w(x » , t) — velocity of point at x^. in + z direction (i. e. , the integral 
of normal "acceleration at the point x A ). "^Solving for a (t) gives 


a (t) = sin 



IT (x A * " w(x A ' 0 ' (x A " X CG } 6 (t) 


( 24 ) 


Using assumption 3, equation (24) becomes 

/ 

*N 


a (t) = sin 


-1 


1 


U. 


I © *i<‘> - w< x A -‘> - <*A - *CG> • <t, l <25) 


By assumption 1, x A - x ^q is known and by assumption 2,^z^ x ^ * s ^ cnown 

for i = 1, 2, .... N. Thus since w(x A ,t) can be measured and 0(t), z(t), 

. . . z N (t) can be found from (22), the right side of ^25) is known when U is 
known. In most cases ot (t) is small so that sin a - a . Then (25) simplifies 
to 


r n 


a (t) = 


U 


I 


i= 1 


X A 


z.(t) - w(x A> t) - (x A - x qq) 0 (t) 


(26) 


The state variables necessary for the equations of motion area; , 0, z^, 
z N , z i , .... z N . Hence if (N + 1) pitch rate sensors are used, 6, z\, 

. . . , zjt can be obtained from (22). Then z\, . . . , zj^ can be obtained by 
integration. The variable a may be computed from (25) lor (26) when applic- 
able I if w(x A ,t) is obtained from a sensor located at x A . 
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The transformation may be written compactly as 



Standard Form for Pitch Equations for an Aeroelastic Vehicle 
The equations of motion can be written as 



( 28 ) 
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15 


and 
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and 



N + 2 columns 


1 o 

i 

V 

K ■ 

0 


• 

« 

1_ 

> 


-F 

T 

-E 




_X_ 



J 


N columns 


N + 2''columns 



rows 


N + 2 rows 


To add the 6 equations to the system (29), let x 


Then 


e 

X 1 


X 2N + 2 


+ 1 0's 


Hx = Kx + 


"6e 

U 


(-6e) 


(30) 


■*N 


6e 


with 


4+ 

H 


(N+3)^ column 


~ 1 0 0 . . 

0 


0 0 . 

• X 

0 1 0 

. . . o 

0 



0 



0 

H 


and K = 

0 

K 


0 



0 


- 
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To put the equations in terms of the sensed variables let 


y 



d G . i = 1. 2, . . . , N+l 

i 

. i = N+2, N+3, .... 2N+2 

i-N-1 

w(x a ,T) 


Then if sin a = a, using (21) and (27), 
x = Ly 

where 



Finally from (30) and (31), 



(31) 


(32) 
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or 


0 


= L 1 H' 1 k L y + L" 1 H' 1 


U 


(-6e) 


(33) 



BLENDING TWO HATE GYRO OUTPUTS 


Two rate gyros are placed on the vehicle in such a way that one is located 
forward and the other aft of the first bending mode antinode. The output of the 
two gyros is automatically blended. The blender is designed so that it will not 
affect the rigid body output but may be chosen to give a positive, negative, or 
zero pickup for the first mode, which ever is desired. 


The block diagram of the blender and logic circuit involved is shown in 
figure 3. The logic circuit consists of two band-pass filters with the pass 
band centered at the first bending frequency, two absolute value units, and a 
servo or electronic multiplier. The output of the blender is denoted by G 0 . 
The outputs of the forward and aft gyros are denoted by Gp and G^, respec- 
tively. Then 


g F " 6 


N 

.©'xo 

1 = 1 


• t <S)' ... V G A 


N 

4 - Y. 


i= 1 


X G, 


z . 

l 


(34) 


From figure 3, it is clear that 

G o ■ KG f + (1 - K)G A = 4 


N 


i= 1 


where 


\. = K (z\ x + 

1 vy x G p 


d - K )^) 1 


k.z. 

i l 


i - 1, 2, .... N. 


(35) 


X G/ 
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Two band-pass filters are used to discriminate all other frequencies from 
that of the£irst mode. In other words, the blender logic is primarily sensi- 
tive to K (^l) 1 and (1 - Zj. The outputs of the band -pass 

filters are then rectified and compared. The difference drives a servo or 
electronic multiplier which will change K until the difference becomes zero, 
i. e. \j - 0. If it is desired to have X.- small but always of one sign, then 

an attenuator could be used as shown in figure 3. 


CONCLUSIONS 

It is shown that, if a finite number of flexure modes are sufficient to 
represent the aeroelastic distortions of a flexible vehicle, the state of the 
vehicle can be measured with commonly available instruments. The method 
requires a number of sensors of the same order of magnitude as the number 
of flexure modes necessary to sufficiently described the aeroelasticity of 
the vehicle. Also it is shown that, by blending the signals from two rate 
gyros, the effects of the first flexure mode can be eliminated. 
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Figure 1. - The flexible vehicle 




Figure 3. - Block diagram of blender and logic circuit 
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